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Abstract 


A deeper understanding of environmental vulnerability is essential for crafting effective sustainability 
strategies. However, data deficiencies and methodological uncertainties often obscure a 
comprehensive grasp of the underlying issues. This study employs Google Earth Engine (GEE)-coded 
satellite imagery and the Getis spatial autocorrelation test to examine spatial patterns of environmental 
vulnerability in 397 Iranian counties from 2011 to 2021, framed within the theoretical lens of 
vulnerability. 

The results of the Getis test reveal heterogeneous distributions and diverse vulnerability patterns across 
Iranian counties in various dimensions, including air pollution, greenhouse gas emissions, changes in 
temperature and precipitation, groundwater depletion, biodiversity threats, and soil erosion, 
particularly in the central, southern, eastern, and northern regions. Emphasizing the role of distinct 
spatial patterns in terms of "sensitivity" and "exposure," this research suggests that environmental 
policies in Iran should be tailored to regional spatial disparities and socio-ecological characteristics. 
This approach can enhance adaptive capacity and mitigate environmental vulnerability in critical 
areas. The use of algorithmic techniques for environmental vulnerability indexing, combined with 
spatial data and population analysis, provides a robust framework for future research. This framework 
can be applied in national and regional planning to alleviate environmental pressures and improve 
sustainability. 

Keywords: Environmental Vulnerability, Google erarth Engine(GEE) ,Spatial Autocorrolation, 
Sustainability Strategies, Regional Disparities, Adaptive Capacity. 
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الگوی فضایی آسیب‌پذیری محبط‌زیست در ایران 
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حکیده 


شناخت سطح خرد آسیب‌پذیری محیط‌زیست از لحاظ فنی و نظری به‌عنوان نقشه راه دقیق در پایداری محیط‌زیست 
عمل می‌کند. باوجوداین ضعف داده و عدم قطعیت روشی فهم درست مسائل ناشی از آن را مغفول گذاشته است. مقاله 
پیش‌رو با استفاده از تصاویر ماهواره‌ای کدنویسی شده گوگل‌ارث‌انجین (GEE)‏ و روش خودهمبستگی فضایی گتیس 
در قالب نظریه‌ی «آسیب‌پذیری» به بررسی و تحلیل الگوهای فضایی آسیب‌پذیری محیط‌زیست ۳۹۷ شهرستان ON pl‏ 
طی‌سال‌های ۱۳۹۰ تا ۱۴۰۰ می‌پردازد. نتایج آزمون گتیس OUS‏ دهنده توزیع نابرابر و الگوهای متنوع آسیب‌پذیری 
شهرستان‌های ایران در ابعاد مختلف ازجمله آلودگی dye‏ گازهای گلخانه‌ای. تغییرات دما و بارش کاهش آب‌های 
زیرزمینی» تهدید تنوع زیستی و فرسایش خاک به‌ویژه در مناطق مرکزی» جنوب. شرق و نوار شمالی خواهد بود. به 
همین ترتیب این پژوهش با ASE‏ بر نقش sla SI‏ فضایی متفاوت در ابعاد «حساسیت» و «قرار گرفتن در معرض 
خطر». پيشنهاد می کند که سیاست‌گذاری‌های محیط‌زیستی در ایران باید با توجه به نابرابری‌های فضایی و متناسب با 
ویژگی‌های جمعیتی و بوم‌شناختی هر منطقه تنظیم شوند. این رویکرد می‌تواند به بهبود ظرفیت سازگاری و کاهش 
آسیب‌پذیری محیط‌زیستی در مناطق بحرانی کمک کند. استفاده از تکنیک‌های الگوریتمی برای شاخص‌سازی 
آسیب‌پذیری محیط‌زیستی و تلفیق داده‌های فضایی با تحلیل‌های جمعیتی. چارچوبی موثر برای پژوهش‌های آینده ارائه 
می‌دهد که می‌تواند در برنامه‌ریزی‌های ملی و منطقه‌ای به‌منظور کاهش فشارهای محیطی و بهبود پایداری محیط‌زیستی 


به کار گرفته شود. 
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صادقی و همکاران» الگوی فضایی آسیب‌پذیری محیط‌زیست در ایران Y‏ 


کلیدواژه‌ها: آسیب‌پذیری محیط‌زیستی ۰ موتور گوگل‌ارث (GEE)‏ « خودهمبستگی مکانی,راهبردهای پایداری, 
تفاوت‌های منطقه‌ای. ظرفیت ساز گاری. 


۱- مقدمه 


یکی از مهم ترین مسائلی که در حال حاضر به‌عنوان یک تهدید جدی برای جوامع مطرح می‌باشد. ریسک‌های 
محیط زیستی ' است که نه‌تنها آرامش و امنیت را از زندگی انسان‌ها دور می‌کند. بلکه موجودیت بشر را نیز 
تهدید می‌کند )2022 Sadeghi, Azizi & Sadeghi,‏ بنابراین» با افزایش آگاهی نسبت به قریب‌الوقوع بودن 
اثرات تغییرات محیط‌زیست بر جمعیت‌های انسانی» برخی از سیاست‌گذاران و اکثر جامعه علمی را برانگیخت 
تا خواهان اقدام موثرتری باشند )2013 (Guzmán & de México,‏ در همین راستا ارزیابی اکوسیستم هزاره 
سازمان ملل" (۲۰۰۵) هشدار داد که «فعالیت‌های انسانی چنان بر عملکرد طبیعی زمین فشار وارد می کند که 
توانایی اکوسیستم‌های سیاره برای be‏ نسل‌های adest‏ را دیگر نمی‌توان بدیهی تلقی کرد» ) Lee & Diop,‏ 
9 چراکه تقاضای انسان در سیاره زمین از عرضه منابع پیشی گرفته و فشار بیش از توان ظرفیت بر 
اکوسیستم‌ها؛ کاهش تنوع زیستی در سطح جهانی و همچنین افزايش سطح منواکسید کربن " در اتمسفر که خود 
باعث افزایش دمای جهانی» تغییرات آب‌وهوا و اسیدی شدن اقیانوس‌ها می‌شود. رخ خواهد داد ) & Guzmán‏ 
.(de México, 2013‏ 

علاوه بر این» همان‌طور که در گزارش معتبر منتشرشده توسط مرکز تاب‌آوری استکهلم " منعکس شد 
پویایی جمعیت‌ها و استفاده نادرست از منابع مادی» انرژی و زیستی زمین. توسط سیستم اقتصادی جهانی در 
حال حاضر از مرزهای سیاره‌ای در سه حوزه (تغییر آب‌وهواء تنوع زیستی و ورودی نیتروژن به بیوسفر) فراتر 
ab)‏ و حداقل در شش حوزه So‏ شناخته‌شده نیز فراتر خواهد رفت )2009 (Rockstróm et al.,‏ در همین 
راستا اکوسیستم‌های طبیعی آن‌چنان به‌وسیله تخریب محیط‌زیست تحت SU‏ قرارگرفته که قابلیت لازم برای 
خود احیائی و سازگاری در آن‌ها به پایین‌ترین حد ممکن رسیده است(2008 (Reed,‏ بنابراین با توجه به موارد 
فوق تخریب محیطزیست به‌عنوان یک مسئله مهم و Sle‏ نگرانی‌های جهانی بسیاری را برانگیخته است. 
بااین وجود به دلیل ضعف داده‌های ثبتی و عدم قطعیت در روش‌های سنجش آسیب‌پذیری امکان ارزیابی 
یکپارچه محیط‌زیست از رهیافت نظری را دشوار کرده است. 


1. Environmental risks 

2. Millennium Ecosystem Assessment 
3. CO2 

4. Stockholm Resilience Center 


۴ نشریه جغرافیا و مخاطرات محیطی» جلد ۰۱۳ شماره ۳ پائیز ۱۳۰۳ 


در همین راستا و از رهیافت نظری محیطزیست می‌توان اشاره کرد یکی از مهم‌ترین دیدگاه‌ها در 
موردحفاظت از محیط‌زیست, نظریه آسیب‌پذیری است؛ چراکه امکان بررسی و ارزیابی آسیب‌پذیری محیط - 
زیست در ابعاد حساسیت " و قرار گرفتن در معرض خطر " را به ورطه تحلیل می‌کشاند. اصطلاح 
«(آسیب‌پذیری ) از کلمه لاتين vulnus‏ به معنای «زخم) گرفته شده است. ازاین‌رو آسیب‌پذیری به معنای واقعی 
Al‏ «توانا یی u^)‏ شدن» است )2006 .(Torabi, Saravani, Sadeghi & Jahanbazian, 2024; Gallopin,‏ سا پر 
تعاریف عمومی این اصطلاح مشابه هستند. به‌عنوان‌مثال de‏ از دست دادن»(1996 L (Cutter,‏ «مستعد بودن 
برای آسیب دیدن»(2006 ,۸0۵0۲ C‏ يا «درجه‌ای که یک سیستم احتمالاً به دلیل در معرض خطر بودن» آسیب 
می‌بیند» )2015 (Turner,‏ آسیب‌پذیری حتی قبل از اینکه کاربرد گسترده‌ای در زمینه‌های بهداشت عمومی. 
توسعه پایدار و تغییرات آب و هوایی پیدا کند. یکی از مفاهیم اساسی در رشته‌های تحقیقی مرتبط با مخاطرات 
طبیعی بوده است )2007 (Turner, 2015 Adger, 2000: Fiissel,‏ 
در ادبیات گسترده آکادمیک» آسیب‌پذیری اصطلاحی است که با انواع پیامدهای خاص مدیریت سوانح و 
به‌طورکلی با خطرات طبیعی مر تبط است )2001 .(Alwang, Siegel & Jorgensen,‏ علاوه بر این» هم جغرافیای 
انسانی و هم بوم‌شناسی انسانی آسیب‌پذیری را به تغییرات محیطی مرتبط می‌داند )2006 (Adger,‏ این در حالی 
است که تحقیقات دیگر آسیب‌پذیری را برحسب سطح خطر» ظرفیت بازیابی و پاسخ به آن تعریف می‌کند؛ 
بنابراین» آسیب‌پذیری علاوه بر معیارهای خطر که با جنبه‌های فیزیکی, اجتماعی و اقتصادی در ارتباط است؛ 
توانایی مقابله با خطرات و شوک‌های مختلف را نیز توصیف می کند )2014 Proag,‏ :1989 ,وتعطه))؛ ¿zl pt‏ 
دو جنبه اساسی در بحث آسیب‌پذیری وجود دارد: نخست. جنبه خارجی که به عناصر ساختاری شامل میزان 
حساسیت و خطر مواجهه )2001 (McCarthy,‏ و دوم جنبه داخلی که به توانایی پاسخگویی و مقابله با عوامل 
تنش‌زا و اقدامات لازم برای غلبه بر آن‌ها اشاره دارد )1989 -Bohle, 2001: Chambers,‏ 
همان‌طور که (Adger, 2006) Sal‏ اشاره دارد. قرار گرفتن در معرض else‏ تنش‌زا. حساسیت و ظرفیت 
سازشی» سه عنصر حیاتی آسیب‌پذیری هستند. سگنوئن و همکاران )2021 (Segnoen et al,‏ نیز ضمن Ast‏ بر 
ویژگی‌های فوق اشاره می‌کنند که آسیب‌پذیری توسط عوامل یا فرآیندهای اجتماعی» اقتصادی» فیزیکی و 
محیطی تعیین می‌شود که حساسیت یک سیستم را نسبت به تأثیر مخاطرات افزایش می‌دهد. قرارگرفتن در 
معرض خطر و حساسیت زمانی که با هم مشاهده می‌شوند. تصویری از تأثیرات بالقوه محیطزیست بر یک 
رو کی فان اچاد کی کل این ارات oss Loud Das aic ve s clap ea Sl‏ را در ظر تی گرا 
بنابراین Lb‏ به‌عنوان مجزا از رتبه‌بندی آسیب‌پذیری کلی در نظر گرفته شوند. هنگامی‌که این سه مؤلفه با هم 
Sensitivity‏ .1 


2. Exposure 
3. Vulnerability 


صادقی و همکاران. الگوی فضایی آسیب‌پذیری محیط‌زیست در ایران ۵ 


در نظر گرفته می‌شوند. می توانند برای ایجاد یک برآورد کلی از آسیب‌پذیری استفاده شوند. در واقع» این مفهوم 
ترکیبی آسیب‌پذیری, ابزار تحلیلی قدرتمندی را برای بررسی و سنجش جدی‌تر آسیب‌پذیری فراهم می‌کند. 
توصیف حالت‌های حساسیت محیط‌زیست به آسیب. شکنند گی» خطر و ظرفیت سازشی در سیستم‌های مختلف 
از سوی دولت‌ها و تصمیم گیرندگان بسیار موردتوجه قرار گرفته است ) Bardsley & Wiseman, 2012 Vogel,‏ 
(Susanne, Roger & Geoffrey, 2012.,‏ 

با توجه به موارد فوق» واکاوی فضایی و زمانی محیط‌زیست و درنتیجه تأیید نظریه آسیب‌پذیری ارائه شده 
توسط col‏ بین دولتی تغییرات اقلیمی ' (۲۰۱۳) در مورد محیط‌زیست ایران می‌تواند به بسیاری از سژالات 
پاسخ دهد و همچنین توجه سیاست گذاران و سرمایه‌گذاران را به وضعیت کنونی محیط زیست جلب کند. علاوه 
بر این» توانایی برنامه‌ریزی و تصمیم‌گیری در مورد مسائل جامعه مثل محیط‌زیست را به سیاست مدران ارائه 
می‌کند. به همین ترتیب» تمرکز اصلی این تحقیق بر نظریه آسیب‌پذیری و جنبه داخلی آن یعنی آسیب‌پذیری 
در ابعاد حساسیت و در معرض خطر بودن و همچنین نقش پویایی جمعیت‌ها در افزایش ظرفیت سازشی 
خواهد بود. ازاین‌رو تحقیق حاضر درصدد است به این سژال که در ابتدا وضعیت محیط زیست ایران در دو 
بعد حساسیت و در معرض خطر بودن چگونه است؟ و UT‏ الگوی فضایی خاصی در این دو بعد آسیب‌پذیری 


وجود دارد؟ پاسخ دهد. 


Jing dip 

۱-۲ موقعیت منطقه موردمطالعه 

منطقه موردمطالعه Gilles‏ شکل ۱ و بر مبنای تقسیمات کشوری در سال ۱۳۹۰ تعداد ۳۹۷ شهرستان ایران در 
طی سال‌های ۱۳۹۰ تا ۱۴۰۰ خواهد بود که بین طول جغرافیایی YO‏ درجه تا ۶۳ درجه شرقی و عرض 
eos P 16 atas.‏ قیال واھ کد ا عفر Cpl‏ کن جیار di‏ ایی عا ها زیت ور 
اقتصادی-اجتماعی مساحتی معادل ۱/۶۲۰ میلیون کیلومترمربع دارد. اقلیم غالب ایران فراخشک و Sis‏ 
بیابانی " (۶۳ درصد کشور) بوده و مقدار متوسط بارند گی آن ۰ میلی‌متر است ) Mohammadi, Karimzadeh,‏ 


.(Pourmanafi & Soltani, 2018 


1 . Intergovernmental Panel on Climate Change 


. مطابق طبقه بندی دومارتون (De Martonne)‏ 
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شکل ۱- موقعیت منطقه موردمطالعه 

Fig.1. Location of the study area 
روش انجام پژوهش‎ ۲-۲ 
x4 پژوهش حاضر از لحاظ هدف. کاربردی و روش آن تحلیل ثانویه است که سطح تحلیل آن مطابق شکل‎ 
تقسیمات کشوری در سال ۱۳۹۰ تعداد ۳۹۷ شهرستان ایران در طی سال‌های ۱۳۹۰ تا ۱۴۰۰ خواهد‎ che بر‎ 
گازهای‎ cd gh داده آلودگی‎ AS ریاست جمهوری (۱۳۹۶) با تمرکز بر داده‌های محبط زیست شامل ۶ مجمو‎ 
d در همین راستا یکی از آلاینده‌های اصلی هوا ذرات معلق‎ (Sadeghi, 2023) ایران موردتوجه قرار گرفت‎ 
به‌وپژه ذرات با قطر آثرودینامکی کمتر از ۲/۵ میکرون " و ۰ میکرون است. که هنوز علی‌رغم تلاش‌های بسیار‎ 
مانده است. اگرچه اندازه گیری‌های مربوط به ایستگاه-‎ BL به‌عنوان یک مشکل فراگیر‎ wal برای کنترل و کاهش‎ 
اما دورسنجی ماهواره‌ای در بسیاری‎ «Raispour, 2021) مطالعه آلودگی 3 اثرات ناشی از این ذرات ناکافی است‎ 
از مطالعات به دلیل در دسترس بودن داده‌ها در یک شرایط با پوشش کامل و یکنواخت جهت دستیابی به توزیع‎ 
محصول متداول‎ 4$ «Han et al., 2015; Kloog et al., 2015: Chudnovsky, 2021 ) بزرگ‌مقیاس ذرات استفاده شد‎ 


1 PM 
2, PM2.5 
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ماهواره‌ای برای برآورد غلظت زمینی ذرات ۲/۵ میکرون و ۱۰ میکرون» عمق نوری هواویزها از طریق ماهواره 
مادیس با وضوح ۰ متر می‌باشد )2021 (Raispour,‏ 

گازهای گلخانه‌ای مجموعه‌ای از گازهای بخار آب. دی اکسیدنیتروژن دی اکسید کربن و متان را در بر 
می گیرند. امروزه به علت تولید بیش‌ازحد گازهای گلخانه‌ای احتمال می‌رود که زمین از حالت استاندارد گرم تر 
شده و باعث» آب شدن یخ‌های قطبی و ایجاد سیل در نقاطی از جهان افزایش فعالیت‌های زمین‌شناختی پوسته 
ne Joys‏ سونامی‌هاء آتشفشان‌ها و طوفان‌ها بیشتر و اسیدی‌تر شدن OT‏ اقیانوس‌ها (و در نتیجه صدمه به 
موجودات آبزی) شود )2011 cpl ply “(Azmi & Motiee,‏ همان‌طور که اشاره شد دو گاز گلخانه‌ای دی اکسید 
کربن و متان از طریق دو ماهواره سنتینل و گوست به‌عنوان مهم‌ترین گازهای مطرح(2013 (IPCC,‏ موردبررسی 
قرار گرفت. 

تغییرات اقلیمی ذیل دو مؤلفه دمای سطح و بارش مدل‌سازی و مورداستفاده قرار می‌گیرد )2013 IPCC,‏ 
از آنجایی دمای سطح زمین یک متغیر کلیدی در آنالیز و مدل‌سازی بیلان انرژی سطح» رطوبت سطحی» تبخیر 
و تعرق و مطالعات تغییر اقلیم در مقیاس‌های گوناگون است ) Sabziparvar & Nazemosadat, 2016: Lafortezza,‏ 
.(Carrus, Sanesi & Davies, 2009‏ می‌توان alan ga:‏ ماهواره‌ای دور سنجی مادیس با رزلوشن ۱۰۰۰ pe‏ در 
بازه‌های زمانی کوتاه نمایشی سینوپتیکی از دمای سطح زمین را ارائه کرد ) & Rozenstein, Qin, Derimian‏ 
(Karmieli, 2014‏ علاوه براین بارش پا بارندگی مهم‌ترین عاملی است که به‌طور مستقیم در چرخه هیدرولوژی 
دخالت دارد )2021 .(Raispour,‏ همان‌طور که اسکندری و همکاران ) Eskandari, Ahmadloo, Pourghasemi,‏ 
(Ahangaran & Rezapour, 2023‏ اشاره می کنند ميانگین بارندگی در کشور Ol pl‏ به‌طور میانگین ۲۴۰ میلی‌متر 
است که از این مقدار حدود ۷۰ درصد آن مستقیماً تبدیل به بخار می‌شود. علاوه براین» برآورد مقدار بارش از 
طریق تصاویر پنج‌روزه ناسا با کیفیت ۵۰۰۰ متر برای پیش‌بینی و شبیه‌سازی تغییرات اقلیمی و مخاطرات طبیعی 
از قبیل سیل. خشکسالی و رانش زمین از اهمیت بالایی برخوردار است )2023 Eskandari et al.,‏ 

افت سطح سفره‌های آب زیرزمینی aS)‏ دومین منبع آب شیرین در دنیا است) به‌عنوان چالشی جدی در 
بیشتر کشورها ایجاد شده است )2003 wil 26 5 Ske, (Cunningham, Sharp, Hiebert & James,‏ چهره‌ی زمین 
را دگرگون خواهد کرد؛ به‌طوری‌که کیفیت آب‌ها را کاهش و هزینه استحصال آن را افزايش می‌دهد ) Issazadeh‏ 
Argany, 2021; Nie et 9‏ &( در چند دهه اخیر در دشت‌های مختلف کشور ایران نیز افت سطح آب زیر 
زمینی شکل گرفته و استفاده از تکنیک دور سنجی و تصاویر ماهواره گریس در زمینه‌های مختلف محیط زیستی 
به‌ویژه منابع آب زیر زمینی گسترش 3b‏ است )2014 Farokhnia & Morid,‏ 

مجموعه‌ای از مطالعات در ارتباط با شاخص پوشش گیاهی نرمال شده نشان می‌دهد که مناطق با پوشش 


گیاهی بیشتر امکان بیشتری را برای زیست انواع گونه‌های جانوری فراهم می‌کند ( De Jong, Verbesselt,‏ 


۸ نشریه جغرافیا و مخاطرات محیطیء جلد ۰۱۳ شماره ۳ پائیز ۱۳۰۳ 


(Schaepman & De Bruin, 2012‏ به همین منظور از تصاویر روزانه مادیس با وضوح ۰ متر و تفاوت بین 
باندهای قرمز و مادون‌قرمز شاخص پوشش گیاهی نرمال شده به‌منظور بررسی ابعاد محیط زیستی پوشش 
گیاهی» تراکم پوشش جنگل‌ها و بازنمایی تنوع زیستی استفاده شد. 

فرسایش خاک" سطح از دست رفتگی خاک توسط عوامل مختلف محیطی نظیر آب‌وهوا شاک توپوگرافی 
و پوشش گیاهی را نشان می‌دهد ) Chen et al., 2019; Thomas, Joseph & Thrivikramji, 2018; Lee, Chu,‏ 
“Guzman & Botero-Acosta, 2021‏ بااین‌حال. مداخلات انسانی از طریق تغییر کاربری زمین ؛ فعالیت‌های 
کشاورزی و ساختمانی می‌تواند این فرآیند را سرعت بخشد ) Vanacker et al., 2019; Diwediga, Le, Agodzo,‏ 
(Tamene & Wala, 2018‏ ازاین‌رو امروزه فرسایش خاک به‌طور گسترده‌ای توسط مدل راسل محاسبه می‌شود 
Paul, Techen, Robinson & Helming, 2019)‏ 

با توجه به موارد فوق و مطابق جدول Y‏ داده‌های مورداستفاده از تصاویر ماهواره‌ای مادیس " (ذرات معلق» 
cles‏ سطح» تنوع زیستی). ماهواره سنتینل و گوست" «گازهای گلخانه‌ای متان و دی‌اکسید کربن)» تصاویر 
ماهواره‌ای GPM‏ و Trmm‏ و شبیه‌سازی پنج‌روزه ناسا" (بارش)» ماهواره گریس GN)”‏ زیر زمینی)» نقشه 
خاک فائو " نقشه مدل رقومی ارتفاع Ol pl‏ و مدل راسل" (فرسایش (SE‏ از طریق کد نویسی در محیط ابری 
گوگل ارث انجین" به دست آمد. به‌طورکلی در بحث داده‌های مورداستفاده امروزه استفاده از تصاویر ماهواره‌ای 
در تکمیل و ارائه داده‌های محیط زیستی به‌طور خاص در پایش آلودگی هوا ) Han et al, 2015: Kloog et al.,‏ 
Chudnovsky, 2021; Karimian et al., 6‏ ;2015( تغییرات اقلیمی ) & Rozenstein et al., 2014; Weng, Fu‏ 


Gao, 2014; Ashouri et al., 2015; Kirschbaum et al., 2017; Yuan et al., 2018; Tang, Clark, Papalexiou, Ma 
«Yoshida et al., 2011; Saxena and Naik, 2019; Schneising et al., 2019) گازهای گلخانه‌ای‎ »& Hong, 2020 


تنوع زیستی ) 2020 «Rannow & Neubert, 2014; Najafi, Darvishsefat, Fatehi & Attarod,‏ منابع آب زیرزمینی 
Feyisa, Meilby, Fensholt & Proud, 2014; Jiao, Zhang & Wang, 2015; Moshayedi & Jahangir, 2021; )‏ 
(Famiglietti et al., 2011‏ و فرسایش خاک )2021 (Mohammadi et al., 2018; Teng et al., 2016; Lee et al.,‏ 


با توجه به تنوع اقلیم‌ها موردتوجه قرار گرفته است. 


. Soil Erosion 

. Land-use change 

. MODIS 

. Sentinel & GOSAT 

. NASA Pentad Precipitation 
. Greece 

FAO 

. RUSSEL 

. Google Earth Engine 
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حدول -į‏ زمان 3 منبع جمع‌آوری داده‌های مو رداستفاده 


Table 1- Time and data collection source 


داده ماهواره (e^‏ جمع‌اوری im‏ زمان 
Time Data collection Source Satellite Data‏ 
Quality‏ 
ذرات معلق "MODIS/006/MCD19A2_GRANULES" Modis‏ 1000 
select('Optical Depth 047? Particulate Matter‏ 
("COPERNICUS/SSP/OFFL/L3. CH4") ۱‏ 
he ۱ Oke‏ ' 
گازهای گلخانه‌ای select('CH4_column_volume_mixing_ratio_d. Sentinel‏ 
D Methane‏ 
ry air) & Greenhouse‏ 1000 
Gases‏ کربن ("COPERNICUS/SSP/OFFL/L3 CO") Gosat‏ 
select('CO_column_number_density’). Carbon‏ 
— بارش r 5-day‏ 
تغییرات اقلیمی ("UCSB-CHG/CHIRPS/PENTAD") Nasa Ran‏ 5000 
Climate‏ 
دما . ("MODIS/006/MOD11A1")‏ 
select(LST_Day_1km). Modis Temp. Change‏ 1000 
آب‌های زیر زمینی S (NASA/GLDAS/V022/CLSM/G025/DA1D')‏ 
select((GWS_tavg') pe‏ 3 
S avg Groundwater‏ 
"MODIS/006/MOD13Q1" Rar gs‏ - 
KANIKI Sg Modis aic‏ ( 250 & 
ae select( ). Biodiversity‏ 
باران ^ = 
S 5000 ("UCSB-CHG/CHIRPS/PENTAD") 2 ue‏ 
e Nasa Rain‏ 
فا ("MODIS/006/MOD13Q1") db‏ 
select(NDVT). Modis Vegetation on‏ 220 
راسل Cover‏ 
Soil erosion‏ توپوگرافی DEM DEM‏ 20 
Topography by the‏ 
RUSLE‏ حفاظت 
DEM DEM model‏ 30 
Conservation‏ 
خاک FAO‏ 
DSMW Soil‏ 30 
Map Soil‏ 


در ادامه برای نشان دادن الگوهای فضایی go‏ >35 در داده‌های تحقیق از آزمون 5 استفاده AS‏ تحلیل 
لکه‌های داغ آماره Getis- Ord)‏ می‌تواند مناطقی را که میانگین محلی آن‌ها بیشتر از میانگین کل است مشخص 
(Getis & Ord, 1996) us‏ این ابزار نه‌تنها برای AU‏ با رد وجود خودهمبستگی فضایی محلی است» بلکه 
اطلاعات قابل نقشه شدن را برای بصری کردن داده‌ها فراهم می‌کند ) Carroll, Reid & Smith, 2008; Hinman,‏ 


2017( ازاین‌رو به‌منظور بررسی ترکیب حوزه فضایی طبق رابطه ۱ از آماره گتیس استفاده AS‏ 
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ei js Wij Xj — WiX 
EU umm RE ۱ رابطه‎ 
2۶ 1 
مجموع وزنی‎ Wi یک وزن مکانی است و مناطق همسایه [به منطقه 1را تعریف می‌کند.‎ Wij در رابطه فوق‎ 
اثبات‌گرایی" برای اثبات نظریه آسیب‌پذیری با بررسی داده‌های‎ Ly است. در گام آخر از دیدگاه‎ Wij مکانی‎ 


محبط زیست در ابعاد حساسیت و در معرض خطر بودن استفاده شد. 


Y‏ نتایج و بحث 
ابتدا یافته‌ها حول چند محور کلیدی به‌صورت توصیفی از روند و الگوی فضایی پارامترهای آلودگی dya‏ 
گازهای گلخانه‌ای» تغییرات اقلیمی. آب‌های زیرزمینی» تنوع زیستی. فرسایش خاک ارائه. سپس نتایج آزمون 
آماری خودهمبستگی فضایی گتیس در بعد میانگین سالانه هر پارامتر (برای نشان دادن شهرستان‌های در معرض 
آسیب‌پذیری و شناسایی محیط‌های آسیب‌پذیر)» عرضه و در انتها میزان رشد سالانه هر پارامتر (برای OLS‏ 


دادن حساسیت) ارائه می‌شود. 


۱-۳ یافته‌های وضعیت محیط زیست ایران 

گردوغبار: 

نتایج شکل OUS Y‏ می‌دهد که به‌طور خاص ذرات Glee‏ ایران در نوار مرزی غرب تا جنوب غرب. جنوب و 
جنوب شرق را پوشش می‌دهد که دریاچه ارومیه و دو کانون مرکزی ذرات Glee‏ در مرکز کشور نیز قابل 
مشاهده است. علاوه براین» نتایج میانگین کشوری ذرات معلق OLS‏ می‌دهد در سال ۱۳۹۰ میانگین سالانه 
شاخحص ذرات Glas‏ ۸۵/۷۲ بود که در سال WAY‏ به مقدار ۷۷/۷۸ و در سال ۱۴۰۰ به مقدار ۷۵/۵۹ کاهش 
می‌یابد. همچنین در سال ۱۳۹۰ شهرستان ارسنجان فارس با میانگین ۱۴۵ میکرون آلوده‌ترین شهرستان ایران و 
شهرستان بهار استان فارس با مقدار میانگین ۱۴ پاک‌ترین شهرستان کشور بوده است. بااین‌حال در سال ۱۴۰۰ 
شهرستان بندرلنگه» استان هرمزگان با مقدار ۱۵۵ آلوده‌ترین و شهرستان بردسیر استان کرمان با میانگین ۲۶/۵۵ 


میکرون پاک‌ترین شهرستان ایران بوده است. 


1. Post positivism 
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۱۴۰۰ شکل ۲- الگوی فضایی و روند تغییرات میانگین ذرات معلق کشور ایران طی سال‌های ۱۳۹۰ تا‎ 
Fig.2. The Spatial pattern and trend of changes in the average particulate matter across 
Iran from 2011 to 2021 


گازهای گلخانه‌ای: 

توزیع گاز گلخانه‌ای دی‌اکسید کربن در کشور Gilles‏ شکل ۳ است. به‌طوری‌که شهرستان‌های نوار غرب تا 
جنوب و جنوب شرق و نوار شمالی کشور به همراه نقاط مرکزی کانون اصلی انتشار گاز گلخانه‌ای کربن پایه 
بوده است. علاوه براین روند سالانه انتشار دی‌اکسید کربن نیز نشان می‌دهد که انتشار آن از سال ۱۳۹۰ با 
میانگین ۱۲/۳۴ به میانگین ۱۳/۲۴ در سال ۱۴۰۰ رسیده است. به‌طورکلی نتایج OUS‏ می‌دهد میانگین انتشار 
دی‌اکسید کربن با گذشت زمان در حال افزايش است و شهرستان‌های نوار شمالی» نوار غرب تا جنوب شرق 
و مرکزی بالاترین میانگین انتشار دی‌اکسید کربن را به خود احتصاص داده و به‌عنوان کانون اصلی انتشار 
دی‌اکسید کربن حائز اهمیت هستند. علاوه براین نتایج انتشار گاز متان (شکل ۳) نیز نشان می‌دهد شهرستان‌های 
مرکزی به سمت شمال شرق و شرق. نوار جنوب غرب تا جنوب و جنوب شرق به همراه نوار شرق تا شمال 
شرق کشور بیشترین مقدار انتشار گاز متان را به خود اختصاص داده‌اند. در مقابل شهرستان‌های نوار شمالی 
شمال غرب و شهرستان‌های غربی کمترین مقدار انتشار را داشته‌اند. نتایج روند سالانه انتشار دی‌اکسید کربن 
نیز OUS‏ می‌دهد که از سال ۱۳۹۰ با مقدار میانگین ۱۹۰۵ تا سال ۱۳۹۴ با مقدار ۱۹۴۱ روند LES!‏ افزایشی و 
بعد از آن تا سال ۱۳۹۷ (تقریباً مصادف با همه گیری کووید 14( با مقدار ۱۸۲۷ انتشار دی‌اکسید کربن به‌صورت 


ميانگین کاهشی و مجدداً روند افزایش می‌یابد؛ به‌طوری که در سال ۱۴۰۰ به مقدار میانگین ۱۹۰۷ می‌رسد. 
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The 4-year average map of CO2 
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شکل ۳- الگوی فضایی و روند تغییرات میانگین گازهای گلخانه‌ای (کربن و متان) کشور ایران طی سال‌های 
۰ تا ۱۴۰۰ 


Fig.3. The spatial pattern and trend of changes in the average greenhouse gases(CO2 & CHa) 
across Iran from 2011 to 2021. 


تغییرات اقلیمی: 

نتایج شکل ۴ OUS‏ می‌دهد به‌طورکلی شهرستان‌های نوار غرب» جنوب غرب. جنوب. جنوب شرق و 
شهرستان‌های مرکزی تا شرق از دمای بیشتری برخوردار هستند. در مقابل. شهرستان‌های شمالی. شمال غرب 
از دمای کمتری برخوردار هستند. نتایج همچنین OUS‏ می‌دهد که میانگین دمای سطح کشور از ۳۱/۳۹ درجه 
سانتی‌گراد در سال ۱۳۹۰ به ۳۳/۸۰ سانتی‌گراد در سال ۱۳۹۶ که گرم‌ترین سال Lp!‏ بوده و ۳۳/۳۹ درجه در 


سال ۰ افزایش awl asl‏ به‌طوری که میانگین دمای کشور ۲ درجه سانتی گرد در طول دوره موردبررسی 


۱۳ 


صادقی و همکاران, الگوی فضایی آسیب‌پذیری محیط‌زیست در Ol pl‏ 


افزایش يافته است. از مهم‌ترین پیامد اصلی افزايش دما در یک منطقه جغرافیایی می‌توان به تغییرات اقلیمی 
گسترده اشاره کرد. علاوه بر این» نتایج شکل ۴ نشان می‌دهد الگوی بارش کشور olal‏ پراکنده و از نابرابری 
گسترده‌ای برخوردار است. به‌طوری‌که در شهرستان‌های نوار شمالی و نوارغربی کشور شاهد میانگین بارش 
خوب و در شهرستان‌های شرقی. جنوب شرق. شمال شرق و مرکز میانگین بارش کمتر از استانداردهای ممکن 
است. علاوه براین روند تغییرات میانگین بارش سالانه نیز OUS‏ می‌دهد که میانگین بارش کل کشور از سال 
۰ با مقدار ۲۳۲ میلی متر به ۳۱۲ میلی متر در سال ۱۳۹۷ افزایش aah‏ و در سال ۱۴۰۰ با مقدار میانگین 


بارش ۱۹۵ میلی متری به پایین تریین حد خود در دوره موردمطالعه رسیده است. 
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The 10-year average surface temperature map 


350 


ap 


ME M E 
mm & 
a 
[5 
ig 
5 E 0 


201 1 201 2201 32014201 5201 6201 7201 8201 920202021 
Year 


gia 


la 


Annua 


روند میانگین ده‌ساله بارش 
The 10-year trend of average precipitation‏ 
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The 10-year trend of average surface temperature 
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نقشه میانگین ده‌ساله بارش 
The 10-year average precipitation map‏ 


شکل ۴- الگو فضایی و روند میانگین تغییرات اقلیمی (دما و بارش) کشور ایران طی سال‌های ۱۳۹۰ تا 


Fig.4. The spatial pattern and trend of average climatic changes (Temperature and Precipitation) 
across Iran from 2011 to 2021 


۱۳۰۳ جغرافیا و مخاطرات محیطیء جلد ۰۱۳ شماره ۲ پائیز‎ 4 pia TZ 


منابع آب زیر زمینی: 

شکل ۵ نشان می‌دهد که در سال ۱۳۹۰ شاخحص آب‌های زیر زمینی به نسبت سال ۱۳۸۹ یعنی سال مبنا ۱۴۵۶ 
میلی متر مثبت بوده و بعد از آن به ۲۳۶۹ میلی متر منفی در سال ۱۳۹۶ می‌رسد. به سبب سیلاب سال ۱۳۹۸ 
شاخحص به ۴۴۲۵ میلی متر افزايش و بعد از آن به ۲۵۴۰ میلی متر منفی در سال ۱۴۰۰ کاهش یافت. quU‏ 
مقایسه شهرستانی نیز نشان می‌دهد در سال ۱۳۹۰ عباس آباد استان مازندران با ۳۳ میلی متر منفی بیشترین 
کاهش و شهرستان لنگرود استان گیلان با ۴۷ میلی متر مثبت بیشترین افزایش را داشته‌اند؛ اما با بررسی این 
نتایج در سال ۱۴۰۰ نتایج می‌توان روند کاهشی آب‌های زیر زمینی را بهتر درک کرد. به‌طوری‌که شهرستان 
تنکابن استان مازندران با ۲۵۵ میلی متر منفی بیشترین کاهش و شهرستان هویزه با ۲۰ میلی متر مثبت بیشترین 
افزایش آب‌های زیر زمینی را داشته که مقایسه مقادیر بیشینه و کمینه شهرستانی نشان از کاهش گسترده آب‌های 
زیر زمین کشور طی دوره موردمطالعه دارد. علاوه بر این مطالعه الگوی فضایی موجود در این رابطه OUS‏ 
می‌دهد کاهش سطح آب‌های زیر زمینی در مناطق شمالی. شمال غرب. جنوب و Sled‏ شرق کشور به واسطه 
کاهش بارش‌ها و افزایش بیش از حد برداشت آب‌های زیر زمینی به یک چالش جدی در کشور تبدیل شده 
است. چراکه کاهش آب‌های زیر زمینی با کاهش رطوبت خاک افزایش تغییرات اقلیمی. فرو نشست گسترده 
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Fig.5. The spatial pattern and trend of average groundwater across Iran from 2011 to 2021 
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صادقی و همکاران. الگوی فضایی آسیب‌پذیری محبط‌زیست در ایران ۱۵ 


نتایج شکل ۶ OUS‏ می‌دهد که نوار شمالی کشور به واسطه پوشش گیاهی نرمال شده بالاتر امکان زیست 
گونه‌های بیشتر گیاهی و جانواری را به نسبت plu‏ مناطق فراهم می‌کند. علاوه براین از شمال غرب تا 
بخش‌های انتهایی جنگل‌های زاگرس در استان لرستان نیز می‌توان پوشش گیاهی نرمال شده مناسبی را هرچند 
ضعیف تر از نوار شمالی مشاهده کرد که شرایط قابل قبولی را برای زیست انواع موجودات زنده فراهم نموده 
است. pls‏ مناطق کشور به‌شدت از ضعف پوشش گیاهی نرمال شده و تنوع زیستی پایین رنج می‌برند. علاوه 
براین همان‌طور که روند تغییرات میانگین سالانه در شکل ۶ نیز OUS‏ می‌دهد در سال ۱۳۹۰ شاخحص پوشش 
گیاهی نرمال شده از مقدار ۰/۱۷ به مقدار ۰/۲۱ در سال ۱۳۹۸ و مقدار ۰/۱۹ در سال ۱۴۰۰ افزایش يافته است. 
به‌طوری که مقایسه نتایج شهرستانی OLE‏ می‌دهد شهرستان فهرج استان کرمان با مقدار ۰/۰۵۷ ضعیف‌ترین 
تنوع زیستی و شهرستان عباس آباد مازندران با مقدار ۰/۶۱ بهترین تنوع زیستی را به خود اختصاص می‌دهد. 
همچنین در سال ۱۴۰۰ شهرستان فهرج با مقدار ۰/۰۶ هنوز ضعیف‌ترین تنوع زیستی و شهرستان عباس آباد با 
افزایش چند صدم درصدی با مقدار ۰/۶۷ بهترین تنوع زیستی و پوشش گیاهی‌نرمال شده را به خود احتصاص 
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The 10-year trend of biodiversity 
۱۴۰۰ شکل ۶- الگو و روند میانگین تغییرات تنوع زیستی کشور ایران طی سال‌های ۱۳۹۰ تا‎ 
Fig.6. The spatial pattern and trend of average biodiversity across Iran from 2011 to 2021 

فرسایش خاک: 
به‌طورکلی الگوی فضایی فرسایش خاک مطابق نتایج شکل ۷ سه کانون نوار شمالی OS‏ نوار غربی و منطقه 
کویری و جنوب شرق کشور را به‌عنوان کانون‌های اصلی فرسایش خاک شناسایی کرد. علاوه براین نتایج روند 
تغییرات فرسایش خاک نیز نشان می‌دهد در سال ۱۳۹۰ soba‏ میانگین ۲۱ تن بر هکتار بر سال فرسایش خاک 
در dl pl‏ وجود دارد که این مقدار به YA‏ تن بر هکتار بر سال برای سال ۱۳۹۷ افزایش بیدا می کند و در سال 


۰ به ۳ تن بر هکتار بر سال می‌رسد که می تواند پیامدهای اقتصادی» اجتماعی و سیاسی گسترده‌ای در پی 


۱۳۰۳ جغرافیا و مخاطرات محیطی» جلد ۰۱۳ شماره ۲ پائیز‎ 4 pia M 
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هکتار بر سال فرسایش خاک وجود داشته که مجموع خاک از دست رفته در این دوره به‌طور متوسط برای هر‎ 
سال معادل ۲۵۴۲۹۲۰۴۴۰ تن است‎ 
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نقشه میانگین ده‌ساله فرسایش خاک 
روند میانگین ده‌ساله فرسایش خاک . انگین aw‏ 
The 10-year trend of average soil erosion‏ 
The 10-year trend of soil erosion‏ 


۱۴۰۰ شکل ۷- الگو و روند میانگین تغییرات فرسایش خاک کشور ایران طی سال‌های ۱۳۹۰ تا‎ 
Fig.7. The spatial pattern and trend of average soil erosion across Iran from 2011 to 2021 

۲-۳ خود همبستگی فضایی گتیس 
در بحث ذرات معلق به‌عنوان یکی از شاخحص‌های GUS‏ سنجش آلودگی هوای کمتر از ۲/۵ میکرون و شاعص 
گرد و غبار با قطر ۱۰ میکرون. نتایج Gilles‏ شکل ۸ در ارتباط با نقشه میانگین سالانه نشان داد که شهرستان‌های 
نوار جنوب غرب تا جنوب و جنوب شرق و شهرستان‌های مرکزی تمرکز بیشترین انتشار ذرات Glee‏ را دارند. 
از انتشار ذرات TU‏ شناسایی می AS‏ علاوه بر این نتایج ميزان رشد انتشارسالانه نیز نشان می دهد که در 
شهرستان‌های غرب و شمال غرب (رنگ آبی) با میزان رشد منفی بین ۳/۸۳ تا ۲/۱۵ شاهد کاهش ذرات Glas‏ 
۵ تا ۲/۸۳ بیشترین میزان رشد ذرات معلق قابل مشاهده است. به همین ترتیب» آزمون گتیس نیز نشان 
دهنده JSS‏ گیری لکه Plo‏ در نوار شمال تا جنوب به‌عنوان کانون اصلی مناطق حساس به آسیب‌پذیری می‌باشد. 
میانگین] که هم رشد انتشار بیشتر و هم میانگین انتشار بیشتری دارند مناطق با ریسک بالا خواهند بود. 
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۱۴۰۰-۱۳۹۰ شکل ۸- نقشه میانگین. میزان رشد و شاخص خودهمبستگی فضایی ذرات معلق طی سال‌های‎ 


Fig.8. The map of the average, growth rate and spatial autocorrelation index of particulate 
matter from 2011 to 2021 


همان‌طور که نتایج شکل ٩‏ نشان می‌دهد در ارتباط با نتایج میانگین گازگلخانه‌ای متان شهرستان‌های مرکزی 
و نوار جنوب غرب تا جنوب و جنوب شرق کانون اصلی انتشار گاز متان و مناطق آسیب‌پذیر را تشکیل 
می‌دهند که نتایج آزمون گتیس نیز تمرکز لکه Flo‏ در این مناطق به‌عنوان شهرستان‌های در معرض آسیب‌پذیر 
را تأیید می‌کند. همچنین در ارتباط با میزان رشد انتشار متان نیز شهرستان‌های شمالی» شمال شرق و مرکز 
بیشترین میزان تغییرات بین (۰/۰۰۰۳ تا ۰/0۰۱۷) به‌عنوان شهرستان‌های حساس به آسیب‌پذیری را نشان می - 
دهد که نتایج آزمون خود همبستگی فضایی گتیس نیز شکل گیری لکه‌های داغ با بیشترین میزان رشد انتشار 
متان در شهرستان‌های نوار شمالی» مرکزی و شمال شرق کشور به‌عنوان مناطق حساس را تأیید می‌کند. علاوه 
براین» شهرستان‌های مرکزی با رشد و میانگین انتشار بالا به‌عنوان مناطق پر مخاطره در انتشار گاز گلخانه‌ای 
متان شناسایی می گردد. 
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Fig.9. The map of the average, growth rate, and spatial autocorrelation index of methane greenhouse 
gas from 2011 to 2021. 


علاوه بر این در رابطه با گاز دی‌اکسید کربن و مطابق نتایج شکل ۱۰ شهرستان‌های مرکزی به سمت شرق» 
جنوب و جنوب شرق با بیشترین میانگین انتشار دی‌اکسید کربن به‌عنوان مناطق در معرض آسیب‌پذیر شناسایی 
شد که نتایج آزمون گتیس نیز شکل‌گیری لکه‌های داغ با تمرکز بیشترین میانگین انتشار دی‌اکسید کربن را در 
این مناطق تأیید می‌کند. همچنین» نقشه میزان رشد تغییرات OUS‏ می‌دهد که بیشترین ميزان رشد تغییرات در 
شهرستان‌های شرقی» نوار جنوب تا غرب و نوار شمالی با میزان رشد بین ۰/۰۰۵۳ تا ۰/۱۲۱ بیشترین ميزان 
انتشار گاز گلخانه‌ای دی‌اکسید کربن را تجربه کرده و به‌عنوان مناطق حساس آسیب‌پذیری معرفی می‌گردد. به 
همین ترتیب نتایج آزمون خود همبستگی فضایی گتیس نیز وجود لکه‌های داغ در نوار جنوب غرب تا جنوب 
و جنوب شرق و شرق را به‌عنوان کانون‌های ple‏ آسیب‌پذیری ناشی از انتشار گاز گلخانه‌ای دیاکسید 


\4 


صادقی و همکاران» الگوی فضایی آسیب‌پذیری محیط‌زیست در Ol pl‏ 


کربن del‏ می‌کند. علاوه بر این شهرستان‌های شرق و جنوب شرق به دلیل برخورداری از میزان رشد و میانگین 
انتشار بالا Ol pea‏ مناطق پر مخاطره شناسایی می گردد. 
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Fig.10. The map of the average, growth rate, and spatial autocorrelation index of CO» greenhouse 
gas from 2011 to 2021. 


در ارتباط با پارامتر تغییرات اقلیمی که به‌وسیله دو مولفه دمای سطح و بارش عملیاتی می‌شود نیز نتایج شکل 


۱ نشان داد در شهرستان‌های شمالی و شمال غرب با دامنه ۱۶/۸۳ تا ۹ عمترین ميانگین دما و شهرستان- 


های نوار جنوب غرب تا جنوب و جنوب شرق» شرق» شمال شرق و مرکز بیشترین میانگین دما را دارند که 


نتایج آزمون خود همبستگی فضایی گتیس نیز با OLS‏ دادن لکه‌های ¿l>‏ در شهرستان‌های غرب. جنوب غرب. 


جنوب. جنوب شرق» شرق و مرکز به‌عنوان کانون اصلی گرماء در معرض آسیب‌پذیر بودن این مناطق را تأیید 


نشریه جغرافیا و مخاطرات «dao‏ جلد OY‏ شماره ۳ پاثیز ۱۴۰۳ 


می کند. علاوه بر این نقشه ميزان رشد دمای سالانه نشان می‌دهد که شهرستان‌های نوار شمالی؛ شمال عرب و 


غرب تمرکز لکه‌های داغ با بیشترین میزان رشد دما به‌عنوان شهرستان‌های مناطق حساس آسیب‌پذیری خواهند 
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Fig.11. The map of the average, growth rate, and spatial autocorrelation index of climate change in 
temperature from 2011 to 2021. 


غربی کشور بیشترین میزان بارش و در شهرستان‌های نوار مرکزی تا جنوب و جنوب شرق کمترین میانگین 


۳۱ 


صادقی و همکاران, الگوی فضایی آسیب‌پذیری محیط‌زیست در Ol pl‏ 


بارش وجود دارد که نتایج آزمون گتیس نیز آسیب‌پذیری شهرستان‌های مرکز تا جنوب و جنوب شرق را 


به‌عنوان لکه داغ با کمترین بارش و شهرستان‌های در معرض آسیب‌پذیر تأیید می‌کند. در مقابل نقشه میزان 


تمرکز لکه plo‏ و میزان رشد منفی بین ۰/۰۳ تا ۰/۰۸ به‌عنوان شهرستان‌های حساس به آسیب‌پذیری Jeb‏ 


می‌شود. به‌طورکلی می‌توان از شهرستان‌های مرکزی با رشد منفی و میانگین کم بارش به‌عنوان شهرستان‌های 
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Fig.12. The map of the average, growth rate, and spatial autocorrelation index of climate change 
in percipitation from 2011 to 2021. 


pia ۳۲‏ 4 جغرافیا و مخاطرات محیطیء جلد ۰۱۳ شماره ۲ پائیز ۱۳۰۳ 


در رابطه با تغییرات آب‌های زیرزمینی نیز مطابق نتایج شکل ۱۳ شهرستان‌های نوار شمالی. مرکز و جنوب 
غرب بیشترین میانگین آب‌های زیر زمینی و شهرستان‌های شمال غرب. شمال شرق و نوار غرب تا جنوب با 
کمترین میانگین آب‌های زیرزمینی شهرستان‌های در معرض آسیب‌پذیری خواهند بود؛ بنابراین» تمرکز لکه داغ 
با کمترین میانگین آب‌های زیرزمینی نیز آسیب‌پذیر بودن این شهرستان‌ها را تأیید می‌کند. علاوه بر اين» نتایج 
میزان رشد سالانه نیز OUS‏ می‌دهد که شهرستان‌های شمالی و شمال غرب و شمال شرق با میزان تغییرات منفی 
بین ۱/۴۵۳ تا ۰۱/۰۹۱۴ بیشترین کاهش منابع آبی زیرزمینی را تجربه کرده‌اند. به‌طوری‌که نتایج آزمون گتیس نیز 
وجود لکه‌های داغ در شهرستان‌های شمالی و شمال غرب به‌عنوان شهرستان‌های حساس به آسیب‌پذیری را 
تأیید می‌کند. به‌طورکلی شهرستان‌های شمال غرب» Sled‏ شرق و جنوب کشور به‌عنوان شهرستان‌های با 
میانگین کم و رشد منفی به‌عنوان شهرستان‌های با مخاطره بالا شناسایی گردید. 
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Fig.13. The map of the average, growth rate, and spatial autocorrelation index of groundwater 
resources from 2011 to 2021 
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صادقی و همکاران. الگوی فضایی Spice‏ محبط‌زیست در ol ul‏ ۳۳ 


پارامتر تغییر انوهی پوشش گیاهی يا پوشش گیاهی استاندارد شده تنوع گونه‌های گیاهی را در یک منطقه 
خاص به تصویر می‌کشد. در همین راستا مطابق نتایج شکل ۱۴ نقشه میانگین تغییر آنوهی پوشش ALS‏ نشان 
می‌دهد شهرستان‌های مرکز به سمت شرق» جنوب و جنوب شرق کمترین میانگین تغییر انوهی پوشش گیاهی 
را دارند. به همین ترتیب. نتایج آزمون گتیس نیز تمرکز لکه‌های Plo‏ را در این مناطق به‌عنوان شهرستان‌های در 
معرض آسیب‌پذیری تأیید می کند. درمقابل نتایج میزان رشد تغییر انوهی پوشش گیاهی نیز نشان می‌دهد 
شهرستان‌های غربی» جلوب شرقی و مرکزی به سمت جنوب با رشد متفی بین MASA‏ تا ۰/۰۱ تمرکز زشد 
منفی و کاهش تغییر انوهی پوشش گیاهی را به‌عنوان شهرستان‌های حساس به آسیب‌پذیری تأیید می‌کند که با 
توجه به میانگین پایین و میزان رشد منفی بالا تمرکز لکه‌های در معرض آسیب‌پذیری و حساس به آسیب‌پذیری 
را شکل می‌دهند. مضاف بر این از شهرستان‌های جنوب غرب و جنوب متمایل به شرق که میانگین پایین و 


رشد مثبت دارند به‌عنوان شهرستان‌های با عملکرد مناسب یاد می‌شود. 
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Fig.14. The map of the average, growth rate, and spatial autocorrelation index of biodiversity 
from 2011 to 2021 


همان‌طور که نتایج شکل ۱۵ در ارتباط با نقشه میانگین فرسایش خاک نشان می‌دهد شهرستان‌های نوار 


شمالی. شمال غرب و غرب بیشترین میانگین فرسایش خاک را دارند که نتایج گتیس نیز در معرض آسیب‌پذیری 


ناشی از فرسایش خاک در نوار شمال کشور را AG‏ می‌کند. همچنین در ارتباط با رشد فرسایش خاک نیز نوار 


غرب تا جنوب غرب جنوب و جنوب شرق بیشترین میزان رشد مثبت فرسایش خاک را دارند که آزمون 


گتیس نیز حساس به آسیب‌پذیر بودن این شهرستان‌ها را تأیید می‌کند. علاوه براین شهرستان‌های غرب و 


جنوب غرب با میانگین بالا و رشد بالا به‌عنوان شهرستان‌های پر مخاطره شناسایی گردید. 
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Fig.15. The map of the average, growth rate, and spatial autocorrelation index of soil erosion 
from 2011 to 2021 


۴ نتیجه گیری 


صادقی و همکاران. الگوی فضایی آسیب‌پذیری محبط‌زیست در ایران YO‏ 


با توجه به اینکه ایستگاه‌های اندازه‌گیری زمینی پراکنده بوده و دسترسی به آمار آن‌ها معمولاً با تأعیر زمانی 
همراه است. وجود روش‌های دیگری برای برآورد داده‌های محیط‌زیستی ضروری به نظر می‌رسد ) Ashouri‏ 
(et al., 2015: Kirschbaum et al., 2017‏ در این راستا اطلاعات ماهواره‌ای و سنجش‌ازدور می‌تواند جایگزین 
و یا مکمل مناسبی برای ایستگاه‌های زمینی در نظر گرفته شود. استفاده از داده‌های ماهواره‌ای امکان برآورد 
پارامترهای محیط زیستی در محدوده وسیعی از سطح زمین» pe pete‏ برای مناطق صعب ‌العبور را امکان‌پذیر 
ساخته که این موضوع باعث توجه محفقان به داده‌های ماهواره‌ای در سطح گسترده شده است ) Ghajarnia,‏ 
(Liaghat & Arasteh, 2015; Kirschbaum et al., 7‏ 

علاوه براین» روش‌های فضایی به دلیل استفاده از مقادیر محلی و عدم تمرکز بر میانگین‌ها قابلیت‌های 
فراوانی را در نشان دادن الگوهای مکانی سطح خرد و نوسانات واقعی مکانی دارند ) Sadeghi, Mohamadi‏ 
(Masiri & mohammadi, 2022‏ به همین دلیل قدرت و رابطه واقعی متغیرها را بهتر نشان می‌دهند؛ Cpl ple‏ در 
این تحقیق از روش خودهمبستگی فضایی گتیس و نظریه آسیب‌پذیری ub‏ دوره ۱۳۹۰ تا ۱۴۰۰ برای ۳۹۷ 
شهرستان dt pl‏ استفاده شد. 

همان‌طور که «IPCC, 2013 Solomon, 2007;) IPCC‏ اد گر )2006 (Adger,‏ و چمبرز )1989 (Chambers,‏ 
اشاره می‌کنند آسیب‌پذیری از سه عنصر Sle‏ قرار گرفتن در معرض» حساسیت و ظرفیت سازگاری تشکیل 
می‌شود. این تعریف چندبعدی آسیب‌پذیری خود در دو جنبه خارجی» شامل pole‏ ساختاری مرتبط با قرار 
گرفتن در معرض das‏ و حساسیت؛ و جنبه داخلی شامل توانایی پاسخگویی و مقابله با عوامل استرس‌زا 
بررسی می‌شود؛ بنابراین با ASE‏ بر جنبه خارجی آسیب‌پذیری وضعیت محیط‌زیست ایران در دو بعد قرار 
گرفتن در معرض خطر و حساسیت. می توان بسترها و زمینه‌های لازم برای ایجاد ظرفیت سازگاری محیط‌زیست 
[به‌طور حاص با تأکید بر نقش پویایی‌های جمعیتی] که همواره مهم‌ترین هدف شناخت محیط‌زیست بوده را 
افزایش داد. 

در همین راستاء نتایج خودهمبستگی فضایی گتیس در ارتباط با بعد در معرض آسیب‌پذیری L]‏ خطر] نشان 
داد که برای شاخص آلودگی هوا شهرستان‌های نوار جنوب غرب تا جنوب و جنوب شرق و شهرستان‌های 
مرکزی» همسو با نتایج احمدی و درویش )2022 (Ahmadi & Darvish,‏ و شایسته و غریبی ) & Shayesteh‏ 
«(Gharibi, 2022‏ برای گاز گلخانه‌ای متان شهرستان‌های مرکزی و نوار جنوب تا جنوب شرق: همسو با عظمی 
و مطیعی )2011 (Azmi & Motiee,‏ و شایسته و غریبی )2022 «Shayesteh & Gharibi,‏ برای گاز گلخانه‌ای 
دی‌اکسید کربن شهرستان‌های مرکزی به سمت شرق تا جنوب و جنوب شرق همسو با شایسته و غریبی 
(Shayesteh & Gharibi, 2022)‏ و قنادی و همکار ان )2022 (Ghannadi, Shahri & Moradi,‏ تغییر ات اقلیمی 
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(Mohammadi et al., 2018)‏ » رمضانی و همکاران( Ramezani Khojeen, Kheirkhah Zarkesh & Daneshkar‏ 
uis 6‏ ات cult‏ ناقی از بارش شهرستان‌هاق gage S a‏ وب کرت همسو میاق 
و همکاران(2018 (Mohammadi et al.,‏ و ذبیحی و همکاران(2015 Zabihi, Sadeghi & Vafakhah,‏ ؛ آب‌های 
زیرزمینی شهرستان‌های QU‏ غرب تا شمال شرق و نوار غرب تا جنوب. همسو با فروتن و 
همکاران(2014 s (Forootan et al.,‏ فرخ نیا و مرید(2014 «Farokhnia & Morid,‏ تنوع زیستی شهرستان‌های 
مرکز به سمت شرق. جنوب و جنوب شرق. همسو با محمدی و همکاران(2018 (Mohammadi et al.,‏ فرسایش 
خاک شهرستان‌های نوار شمالی» شمال غرب و غرب. همسو با محمدی و همکاران(2018 (Mohammadi et al.,‏ 
با میانگین بالای تخریب در معرض آسیب‌پذیری L]‏ خطر] قرار دارند. 

علاوه quoi‏ نتایج خودهمبستگی فضایی گتیس در ارتباط با بعد حساسیت نیز تأیید می‌کند که برای شاخص 
آلودگی هوا شهرستان‌های جنوب مرکز و شمال به‌صورت یک خط ممتد؛ برای گاز گلخانه‌ای متان شهرستان- 
های شمالی» شمال شرق و مرکز؛ برای گاز گلخانه‌ای دیاکسید کربن شهرستان‌های شرق, نوار جنوب تا غرب 
ووا ما رای کرات ایی اتی ار ا شیرسان‌های تزا یمان ال pee‏ کر برای رات 
اقلیمی ناشی از بارش شهرستان‌های شمالی» شمال غرب مرکز تا شرق؛ برای آب‌های زیرزمینی شهرستان‌های 
شمالی و شمال غرب؛ برای تنوع زیستی شهرستان‌های مرکزی به سمت جنوب؛ برای فرسایش خاک شهرستان- 
های نوار غرب تا جنوب غرب و جنوب شرق با رشد CVE‏ تخریب در وضعیت حساس آسیب‌پذیری قرار 
دارند. 

به‌طورکلی مقایسه نتایج آزمون‌های خودهمبستگی فضایی ابعاد آسیب‌پذیری حساسیت و قرار گرفتن در 
معرض خطر نیز OUS‏ می‌دهد که برای شاخص آلودگی هوا شهرستان‌های جنوبی و مرکزی به سمت شمال؛ 
برای گاز گلخانه‌ای دی‌اکسید کربن شهرستان‌های مرکزی؛ برای تغییرات اقلیمی ناشی از دما شهرستان‌های 
شرق» جنوب شرق؛ برای تغیبرات اقلیمی ناشی از بارش شهرستان‌های مرکزی؛ برای آب‌های زیرزمینی 
شهرستان‌های شمال غرب. شمال شرق و جنوب کشور؛ برای تنوع زیستی شهرستان‌های مرکزی به سمت 
جنوب. جنوب و شرق؛ برای فرسایش خاک شهرستان‌های غرب و جنوب غرب که هم رشد بالا و هم میانگین 
تخریب بالای دارند در دسته شهرستان‌های با آسیب‌پذیری حاد دسته‌بندی می‌شوند؛ بنابراین می‌توان اشاره کرد 
که نتایج تأییدکننده آسیب‌پذیری در ابعاد حساسیت و قرار گرفتن در معرض Jas‏ محیط‌زیست ایران به‌صورت 
الگوهای نابرابر فضایی همسو با مطالعات محمدی و همکاران )2018 «(Mohammadi et al.,‏ شایسته و 
«(Shayesteh & Gharibi, 2022) ga‏ قنادی و همکاران(2022 «(Ghannadi et al.,‏ فرخ نیا و مرید( Farokhnia‏ 


Morid, 2014‏ &( است. 
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با مشخص شدن الگوهای فضایی متفاوت ابعاد حساسیت و در معرض خطر بودن آسیب‌پذیری می‌توان از 
مجموعه پیشنهادهای سیاستی برای کاهش آسیب‌پذیری و افزایش ظرفیت سازشی سخن گفت. به همین ترتیب 
در ذیل پيشنهادها این پژوهش حول چند محور کلیدی دسته‌بندی می‌شود: نخست. نقش سیاست‌های سطح 
خرد با توجه به الگوهای نابرابر شهرستانی جهت سیاست گذاری مرتبط با الگوهای متفاوت هر شهرستان در 
اولویت است. دوم. نقش سیاست‌های دوگانه برای شهرستان‌های در معرض خطر و شهرستان‌های حساس 
حائز اهمیت است. به‌طوری که شهرستان‌های در معرض حساسیت. ASE‏ سياست‌ها بر پیشگیری و در 
شیر cela a‏ در سر خط ر ترک cal‏ بر شیاست‌های: ate bo PET cul cull‏ اک این 
مطالعه بخشی از یک پژوهش گسترده‌تر )2023 (Sadeghi,‏ است با استناد به نتایج پژوهش مذ کور می توان به‌طور 
خلاصه بر رابطه جمعیت و محیط‌زیست به‌عنوان یکی از ملفه‌های مهم افزایش ظرفیت سازشی تأکید کرد. 
جونز و پرستون )2011 Jones & Preston,‏ اشاره می‌کند تعاملات جمعیت و محیط زیست یک میدان مطالعاتی 
به‌طور فزاینده مهم است که نیازمند توجه فزاینده است. هنگام توجه به مسائل محیط‌زیستی اگرچه مسیرهای 
شناخته‌شده شامل نژاد. قومیت. پویایی جنسیتی رشد اقتصادی و همچنین دسترسی به شبکه‌های ایمنی. 
اجتماعی و خدماتی که برای تأمین معیشت ضروری هستند. به‌سرعت دنبال شود. بااین و جود برخی از جنبه‌های 
پوبایی جمعیت. مانند مهاجرت. شهرنشینی. تحولات خانوار تحصیلات و ساختار سنی به‌طور مستقیم با 
ساز کاریق مرتبط هستند )2022 (Sadeghi, 2023; Abbasi Shavazi & Sadeghi,‏ ازاین روء تجزیه‌و تحلیل Sis‏ - 
ها و پویایی جمعیت می‌تواند ابزار قدرتمندی برای برنامه‌ریزی سازگاری و ظرفیت سازشی ایجاد AS‏ 
(Guzmán & de México, 2013)‏ شربین و همکاران )2007 (Sherbinin, Carr, Cassels & Jiang,‏ استدلال 
cl‏ علاقه به موضوع LEU‏ تیبرت tud patho‏ دز gh gud die‏ انو آما el‏ کرک 
پوبایی جمعیت شکاف قابل‌توجهی در توسعه و اجرای پروژه‌های سازگاری ایجاد می‌کند. در همین راستا؛ 
توجه مناسب به رشد. ترکیب و توزیع جمعیت که ساز و کار آن‌ها در پژوهش صادقی Sadeghi,2023)‏ بررسی 
گردید برای درک چگونگی توزیع آسیب‌پذیری محیط‌زیست در بین گروه‌های مختلف مردم حیاتی است (برای 
اطلاعات بیشتر به رساله دکتری 2023 Sadeghi,‏ مراجعه شود). 

در یک نتیجه گیری کلی» محیط‌زیست ایران در پارامترهای آلودگی هواء تغییرات اقلیمی» گازهای گلخانه‌ای, 
منابع آب زیرزمینی» تنوع زیستی و فرسایش خاک از الگوی فضایی متفاوت حساسیت و قرار گرفتن در معرض 
آسیب‌پذیری برخوردار است که آزمون آماری گتیس نیز تأییدکننده الگوهای متفاوت فضایی آسیب‌پذیری 
خواهد oy‏ لذا در بستر توسعه اقتصادی و اجتماعی توجه به ابعاد آسیب‌پذیری و اولویت‌بندی شهرستان‌ها با 
توجه به درجه قرار گرفتن در معرض خطر و حساسیت می‌تواند تأثیر واقعی مؤلفه‌های مؤثر بر ظرفیت سازگاری 
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آسیب‌پذیری» استفاده PE 5l‏ ماهواره‌ای» شاخص سازی ترکیبی ابعاد محبط زیست L‏ استفاده از روش‌های 
الگوریتمی تسلط تصادفی برای تعیین وزن ابعاد گوناگون در سنجش آسیب‌پذیری و ASU‏ بر نقش پویایی 


جمعیت در افزایش ظرفیت سازگاری محبط زیست از مهم‌ترین پیشنهادهای این پژوهش برای مطالعات NUNT‏ 


است. 


تقدیر و تشکر: این مقاله مستخرج از رساله دکتری است که با حمایت دانشکده علوم اجتماعی دانشگاه تهران 


^c‏ شده است. از این دانشکده E‏ دانشگاه تهران p‏ 3 قدردانی می‌شود. 
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